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Abstract 

Thermal instability is one of the most important processes in the formation of 
clumpy substructure in magnetic molecular clouds. On the other hand, ambipolar 
diffusion, or ion-neutral friction, has long been thought to be an important energy 
dissipation mechanism in these clouds. Thus, we would interested to investigate the 
effect of ambipolar diffusion on the thermal instability and formation of clumps in the 
magnetic molecular clouds. For this purpose, in the first step, we turn our attention 
to the linear perturbation stage. In this way, we obtain a non-dimensional character- 
istic equation which reduces to the prior characteristic equation in the absence of the 
magnetic field and ambipolar diffusion. With numerical manipulation of this charac- 
teristic equation, we conclude that there are solutions where the thermal instability 
allows compression along the magnetic field but not perpendicular to it. We infer that 
this aspect might be an evidence in formation of observed disc-like (oblate) clumps in 
magnetic molecular clouds. 
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1 Introduction 

An almost thorough analysis of linear stage of thermal instability was given in a well-known 
paper by Field(1965). He showed that thermal instability can lead to the rapid growth of 
density perturbations from infinitesimal 5p to nonlinear amplitudes on a cooling time-scale, 
t c , in which for typical conditions in the interstellar medium (ISM) is short compared to the 
dynamical time-scale t&. Of the purely thermal modes, the most relevant for the ISM is the 
isobaric mode which has been discussed in terms of the formation of distinct phases of the 
ISM (Field, Goldsmith & Habing 1969) and the formation of protostars in a cooling ISM 
(Schwarz, McCray & Stein 1972). The isentropic mode has also found application in the 
ISM and has been discussed with regard to the amplification of acoustic waves in a warm 
ISM (30 < T < 90 K) (Flannery & Press 1979). A more detailed investigation of the growth 
of condensations in cooling regions has been presented by Balbus(1986,1991) who examined 
the effect of magnetic field. 

High resolution studies of the magnetic molecular clouds, reveal that they have internal 
structures on all scales and are typically clumpy or filamentary (Falgarone, Puget & Perault 
1992, Langer et al 1995), with prolate and oblate (disc-like) clumps (Ryden 1996). Gammie 
et al (2003) have recently studied the three dimensional analogs of clumps. They have con- 
cluded that nearly 90% of the clumps are prolate and 10% of them are oblate. The origin 
and shape of these clumps is a disputable issue. Thermal instability and turbulence may be 
two responsible parameters. 
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In molecular clouds, the dispersion velocity inferred from molecular line width is often larger 
than the gas sound speed inferred from transition temperatures (Solomon et al 1987). MHD 
turbulence may be responsible for the stirring of these clouds (Arons & Max 1975). because 
of these turbulent motions, molecular clouds must be transient structures, and are probably 
dispersed after not much more than ~ 10 7 yr (Larson 1981). Since cooling time-scale of 
molecular clouds is approximately ~ 10 3 — 10 4 yr (Gilden 1984), thermal instability may be 
a coordinated trigger mechanism for clump formation. Turbulence, in the second stage, can 
deform these small-scale clumps in shape and orient them relative to the magnetic fields. 
Gilden(1984) calculated the net cooling function for molecular clouds and found that in 
environments where CO cooling dominates, molecular gas may be thermally unstable. He 
suggested that thermal instability may be an important source of small-scale clumps in 
fully molecular clouds. Burkert & Lin(2000) have recently proposed that clumpyness in 
cold clouds arises naturally from their formation through a cooling instability which acts 
on time-scales that can be much shorter than the dynamical time-scale of the cloud. Af- 
terward, Gomez-Pelaez & Moreno-Insertis(2002) have investigated the effect of self-gravity 
and thermal conduction on a cooling and expanding medium. They classified importance of 
various physical processes including self-gravity, background expansion, cooling, and thermal 
conduction according to their relative time-scales. 

On the other hand, observations establish that magnetic fields play an important role in 
shaping the structure and dynamics of molecular clouds. Especially, ambipolar drift, or ion- 
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neutral friction, has long been thought to be an important energy dissipation mechanism in 
magnetic molecular clouds (Scalo 1977). If A represents the characteristic dimension over 
which the magnetic field varies (the wavelength of perturbation), the time-scale of ambipolar 
diffusion in a typical molecular cloud may be approximated as tad ~ 10 8 ^^pc)^ r (Shu 1991). 
We expect that for a critical wavelength (~ O.Olpc) which tad ~ t c , ambipolar diffusion may 
be important. Thus, in view of this importance, we would interested to investigate the effect 
of it on thermal instability and formation of small-scale clumps in the magnetic molecular 
clouds. We suggest that shape of oblate clumps is formed from their early stage of evolution, 
via thermal instability. 

For this purpose, in the first step, we turn our attention to the linear stage and neglect the 
effect of self-gravity and background expansion or contraction. Section II of this paper devel- 
ops the theory of linear thermal instability in the presence of ambipolar diffusion. In §111, we 
obtain a non-dimensional characteristic equation which can reduce to the Field's character- 
istic equation in the absence of ambipolar diffusion. Then, we discuss about the domains of 
stability (or instability) of this characteristic equation. We find a critical wavelength which 
the effect of ambipolar diffusion is very important and small-scale disc-like clumps can be 
formed. Finally, a conclusion is given in §IV. 
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2 The Linearized Equations 

A molecular cloud gas includes neutral atoms and molecules, atomic and molecular ions, and 
electrons, which are the primary current carriers. Since significant charge separation can not 
be sustained on the time-scale of interest, so the electrons and ions move together. 
In principle, the ion velocity Vi and the neutral velocity v n should be determined by solving 
separate fluid equations for these species (Draine 1986), including their coupling by collision 
processes. But, in our interested time-scale of cooling (10 3 — 10 4 yr, Gilden 1984), two fluids 
of ion and neutral are well coupled by together, thus we can use the basic equations as follows 
(Shu 1991) 

g + „V.„ = (1) 
1 dp 7 P« +(<! _ v . (JfvI , ) = (3) 



7 — 1 dt 7 — 1 p dt 
^ + B(V • v) - (B ■ V)v = V x { x [B x (V x B)]} (4) 

P --pT = (5) 
/i 

where d/dt = d/dt + v • V is the Lagrangian time derivative and K is the coefficient of 
thermal conduction. 7 is the polytropic index of the ideal gas, p is the mean atomic mass 
per particle, R is the universal gas constant, and i] = < ^ n °£ > where v in is the ion-neutral 
relative velocity with impinging cross section Oi n . In writing the above equations, we used 
the relation pi = tp v n between ion density and neutral density (Nakano 1980). Since, the ion 



density is much less than the neutral density, in a good approximation we estimate 

P = Pn + Pi ~ Pn- (6) 

f2(p, T) is the net cooling function (erg. sec' 1 .gr~ v ) as follows 

n(p,T) = A(p,T)-T tot (7) 

where T tot is the total heating rate and A(p, T) is the cooling rate which can be written as 
(Goldsmith & Langer 1978) 

A(p,T) = A p 5 T? (8) 

where A , 5, and (3 are constants. The range of (3 is 1.4 to 2.9. The constant 5 is greater 
than zero for optically thin case and less than zero for optically thick case. Models of molec- 
ular clouds identify several different heating mechanisms; cosmic rays, H 2 formation, H 2 
dissociation, grain photoelectrons, collisions with warm dust, gravitational contraction, and 
ambipolar diffusion. The rate for these processes are largely unknown. Sticking efficiency 
on grains, grain composition and lattice structures, cosmic ray spectra and flux, efficiency of 
cosmic ray penetration into clouds, magnetic field strengths and geometry, and the fractional 
ionization are a few disputable parameters. In this paper we consider the heating rates of 
cosmic rays, H 2 formation, H 2 dissociation, grain photoelectrons, and collisions with warm 
dust as a constant T (Glassgold & Langer 1974, Goldsmith & Langer 1978). Also, we dis- 
regard gravitational heating rate, since we neglect self-gravity and background contraction. 
The heating of the gas by magnetic ion-neutral slip is discussed in detail by Scalo(1977); our 
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simple estimate of this heating rate is as follows 

Y AD = V ep»v 2 d (9) 

where Vd is the drift velocity of the ions relative to the neutral. With above explana- 
tions/approximations and with introduction T' = rjev^, in a good manner we choose the net 
cooling function as follows 

Q(p,T) = A p s T f3 -r -r' oP \ (io) 

In the local homogeneous equilibrium state, we have p = po,T = T ,p = p , B = B , v — 0, 
and Q(po,T ) = 0. We assume perturbations of the form 

A(r,t) = Aiexp(ht + ik-r). (11) 

Then the linearized equations are 

hpi + ip k • v± — (12) 

hp oVl + ik Pl + i(B • Bx) A - i(k • B )5i = (13) 

47T 47T 



i-Pi - 7 ^ + A> n ^i + Aj^tTi + Kk T\ = (14) 

7-1 (7 - l)p 

/iBi + iB (k ■ - i(k ■ B )v 1 = ik x { \— x [B x (ik x Bi)]} (15) 

4nr]epQ 

v- P i-¥ = ° < 16 > 

Po Po Jo 
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where Q p = (dVL/dp)T and VLt = (dVl/dT) p are evaluated for the equilibrium state. 
We introduce the coordinate system u x ,u y , u z specified by 

Bo B x k 

tt, = — - , u v = — — , U x = U v X u z . (17) 

B ' y |B x fc| ' * v * v ) 

Equations (13) and (15) may be used to uncouple V\ y - the perturbed velocity in the plane 
perpendicular to both B and k- from the reminder of the problem. Disturbances perpen- 
dicular to the (B — fc)-plane have a solution of the form 

(k\\a) 2 , r . k\\a , 9l i , . 

where a is the Alfven velocity and k\\ = k ■ u z . Thus, in long wavelength perturbations 
that k\\a < 2r]epQ, waves are damped with ion-neutral friction and in short wavelength 
perturbations that k\\a > 2i]epQ, Alfven waves can not be existed. 

The motion in the other modes are constrained to the x — z-plane, and are governed by the 
characteristic equation, 

h 5 + {c s k T + /i')h 4 + [k 2 {c 2 s + a 2 ) + h'k T c s ]h 3 + {k 2 [c 3 s {k T - k p )h + a 2 c s k T ] + c 2 h'k 2 }h 2 

+ [c 2 s a 2 k 4 cos 2 6 + k 2 c 3 s h'{k T - Ag/ 7 ]h + (k T - k p )c 3 s a 2 k 4 cos 2 0/ 7 = (19) 
where c s is the Laplacian speed of sound, 9 is the angle between k and B , and 

M7-i)PDnp ^ = M7-i)^ + Mi z i)jr A;2 (20) 

Rc s T Rc s Rc s po 

are wavenumbers of sound waves whose angular frequencies are numerically equal to the 
growth rates of isothermal and isochoric perturbations, respectively, h! = k 2 a 2 /r]epQ is the 
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effect of ion-neutral slip. If we neglect the effect of the magnetic field and ambipolar diffusion 
(a = 0, h' — 0), the characteristic equation (19) reduces to the equation (15) of Field(1965). 

3 Domains of Stability 

By introducing the non-dimensional quantities, 

h k kr ,a,, ^ h' 

y =k^ ' a " = ~k ' * T = T ' ' D =k^ (21) 

we can write the characteristic equation in the following form, 

y 5 + (<t t + £>)y 4 + (1 + a + a T D)y 3 + [7 _1 ((Tt - <r p ) + aa T + D]y 2 

+[a cos 2 6 + -f'\a T - a p )D]y + 7^ 1 a(cr T - cr p ) cos 2 9 = (22) 

so that for each # we have four free parameters that consist of ax, cr p , a, and D. We want 
to study the effect of ambipolar diffusion on stable region of the a T — <? p plane. For this 
purpose, we use the Laguerre's method for finding the roots of the characteristic equation. 
First we consider the problem in the absence of the magnetic field (a — 0, D — 0). The stable 
region of this case is shown in Fig. I. This result had been derived by Field(1965). Now, we 
are interested to consider the effect of the magnetic field. For this purpose, we must consider 
the ambipolar diffusion because of small ion density in magnetic molecular clouds. In this 
case, we must break a lance to the complete characteristic equation (22) for different values 
of cry, <7 p , cc, D, and 9. The stable regions of these typical instances are shown in Fig. 2. In 

10 



this case, the line OA in Fig. 1 is unchanged, corresponding to the breaking of the magnetic 
pressure for the reason of ion- neutral slipping and smallness of ions. But, the line OB in 
Fig. 1 is brought down, corresponding to the dissipating ion-neutral slip heating during the 
compression phase of the wave. Thus, ambipolar diffusion can stabilize the medium so that 
it's maximum effect is occurred at 6 — ir/2. 

Inserting the net cooling function, Equ.(lO), into the definitions of a p and <tt, we get 

°>=%i^™- v ® < 23 > 
^=S" A( "»- t » )[1 - ( t' 2 i (24 > 

where £ is the ratio of ambipolar diffusion heating rate to the cooling rate as 
and Aq is a defined wavelength as follows 



Ao^J f To (26) 
V (Jp A(p , T ) 



We separate two cases as follows 



1. 5 > £f which <jp > that is upwards of the a T — a p plane. In this case, the medium is 
optically thin. As shown in Fig. 3(a), two regions of the — a p plane are separated 
by a critical wavelength 
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If the wavelength of perturbation is greater than this critical value, medium is stable. 
If A < A cl , the magnetic molecular cloud is unstable and a spherical clump can be 
formed. 

2. 5 < £f which a p < that is downwards of the a T — a p plane. The optically thick 
molecular clouds set in this case. As shown in the Fig. 3(b), three regions of the 
&t — v P plane are separated by two critical wavelengths 



If the wavelength of perturbation is greater than A c3 , the medium is stable. If A < A c2 , 
the molecular cloud is unstable and a spherical clump can be formed. Depend on the 
values of D and a, we have semi-stable regions between A c2 and A c3 , which the thermal 
instability allows compression along the magnetic field but not perpendicular to it. 

The increased area of semi-stable region (shaded areas of Fig. 2) as a function of D for three 
typical values of a and for 6 = n/2 is plotted in Fig. 4. According to this figure, maximum 
of stability is occurred at a typical D m . Thus, we may define a critical wavenumber 



At small wavelengths which k ^> k c , ambipolar diffusion can break the effect of the magnetic 
field on the whole matter, thus the value of A s is zero. On the other hand, at very large 
wavelengths which k <C k c , ambipolar diffusion time-scale is very greater than the cooling 
time-scale, therefore we can neglect its effect , thus the value of A s must be zero, too. 




(28) 



k c — 



a 2 



D 



(29) 
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4 Conclusion 

We have carried out a systematic analysis of the linear thermal instability of a locally uniform 
magnetic molecular cloud which, in the perturbed state, is undergoing ambipolar diffusion. 
Although thermal instability proceeds faster than dynamical processes such as turbulence, 
its growth rate is determined by the local cooling rate. We choose a simple parametric 
net cooling function and discuss about its different parameters for unstability and clump 
formation. The small perturbation problem yields a complete characteristic equation that 
in the absence of the magnetic field and ambipolar diffusion, reduces to the prior results of 
the linear thermal instability. We have used the Laguerre's method for finding the roots of 
this characteristic equation. 

The stable region by neglecting the magnetic field is shown in Fig. 1, while, Fig. 2 displays 
the stability region for typical values of the magnetic field (a) and the ambipolar diffusion 
strength (D). Comparison of these figures indicate that ambipolar diffusion can stabilize the 
medium, in this manner that its maximum stabilization is occurred perpendicular to the 
magnetic field {6 = 7r/2). Thus, including the magnetic field and considering the ambipolar 
diffusion, divides the <Jt — v P plane in three regions: stable region, unstable region, and 
semi-stable region. 

By inserting the parametric general form of net cooling function into the definitions of ot and 
Op, we find critical wavelengths which divide different cases of stability, instability, and semi- 
stability of the ot — o p plane according to the wavelength of perturbation. If the physical 
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parameters of the molecular cloud (or the wavelength of perturbation) settle on the unstable 
region of Fig. 3, a spherical clump must be formed. On the other hand, if its parameters or 
the wavelength of perturbation settle on the semi-stable region, thermal instability allows 
compression along the local magnetic field but not perpendicular to it. Therefore, including 
the magnetic field and ambipolar diffusion may be an evidence in formation of small-scale 
disc-like clumps in magnetic molecular clouds. 

Since ambipolar diffusion time-scale depends on the wavelength of perturbation, we find a 
critical wavenumber which the effect of ambipolar diffusion for stabilizing the medium is 
very important. 

We have assumed a uniform background. In spite of this, our results are applicable locally in 
a non-uniform background if the perturbation wavelength is much less than the macroscopic 
variation length of the unperturbed quantities. 

We neglected the interaction and merging of the clumps. These processes become important 
for the subsequent evolution. We also neglected the effect of self-gravity and contraction or 
expansion of the background. They will be considered in the subsequent papers. 
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Figure 1: Region of stability in the ot — <J P plane in the absence of the magnetic field 
( a = 0,D = 0). 
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Figure 2: Regions of stability in the ot — c P plane for presence of magnetic field including 
ambipolar diffusion («^0,D^0) for three values of 9. The shaded areas are semi-stable 
regions where a disc-like clump can be formed. 
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Figure 3: Different cases of stability, unstability, and semi-stability of the ot — cr p plane 
in terms of A, wavelength of perturbation. Semi-stable regions of figure (b) depend on the 
magnetic field (a) and the ambipolar diffusion strength (D). 
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Figure 4: The area of semi-stable region, as a function of D for three typical values of a. 
Maximum of semi-stable region is occurred at a typical D m . 
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